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1. Introduction 

The modem theory of strong interactions, QCD (Quantum Chro mo-Dynamics) [|I|, |3|], has a 
simple structure but a very rich dynamical content (for books and reviews, see for example [||, 
It gives rise to a complex spectrum of hadrons, implies the striking properties of confinement and 
asymptotic freedom, is endowed with an approximate chiral symmetry which is spontaneously 
broken, has a highly non trivial topological vacuum structure (instantons, U{\)a symmetry break- 
ing, strong CP violation (which is a problematic item in QCD possibly connected with new physics, 
like axions, ...), an intriguing phase transition diagram (colour deconfinement, quark-gluon plasma, 
chiral symmetry restoration, colour superconductivity, ...)■ 

How do we get testable predictions from QCD? On the one hand there aie non perturbative 
methods. The most important at present is the technique of lattice simulations (for a recent review, 
see ref. [^): it is based on first principles, it has produced very valuable results on confinement, 
phase transitions, bound states, hadronic matrix elements and so on, and it is by now an estab- 
lished basic tool. The main limitation is from computing power and therefore there is continuous 
progress and a lot of good perspectives for the future. Another class of approaches is based on 
effective lagrangians which provide simpler approximations than the full theory, but only valid in 
some definite domain of physical conditions. Typically at energies below a given scale L particles 
with mass larger than L cannot be produced and thus only contribute short distance effects as vir- 
tual states in loops. Under suitable conditions one can write down a simplified effective lagrangian 
where the heavy fields have been eliminated (one says "integrated out"). Virtual heavy particle 
short distance effects are absorbed into the coefficients of the various operators in the effective 
Lagrangian. These coefficients are determined in a matching procedure, by requiring that the ef- 
fective theory reproduces the matrix elements of the full theory up to power corrections. Important 
examples of effective theories are chiral lagrangians, based on soft pion theorems and valid for 
suitable processes at energies below 1 GeV (for a recent, concise review see ref. and references 
therein); heavy quark effective theories [^, obtained from expanding in inverse powers of the heavy 
quark mass, mainly important for the study of b and, to less accuracy, c decays (for reviews, see, for 
example, ref. [|l0[]); soft-coUinear effective theories (SCET) [|ll|], valid for processes where quaifc 
with energy much larger than their mass appeai\ On the other hand, the perturbative approach, 
based on asymptotic freedom, still by far remains the main quantitative connection to experiment, 
due to its wide range of applicability to all sorts of "hard" processes. 

The crucial ingredient of perturbative QCD is the "running" coupling 0:^(2), the scale depen- 
dent effective coupling that decreases logarithmically with increasing values of the energy scale 
Q and tends towards zero at infinity ("asymptotic freedom"). At large scales the detailed form of 
the running is prescribed by the QCD theory. The running coupling is a property of the theory 
and it enters in the perturbative expansion of any hard process. Some particularly simple and clear 
hard processes must be used to measure the running coupling at a suitable energy scale. Once the 
running coupling is measured it can be used to predict all observables related to a variety of hard 
processes. In particular- QCD plays a crucial role in the physics at hadron colliders where many 
important processes, like jets at large pT, heavy quark pair production, Higgs production and many 
more have rates that start at order a^. It is clear that the measurements of a, are very important to 
work out the relevant predictions of the theory. The error on the measured value of a, is a main 
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component of the theoretical error in all perturbative QCD predictions. In this article after recall- 
ing the basic definitions and facts about the running coupling we go to our main focus which is a 
critical discussion of the methods for precisely measuring and the coixesponding results. 

2. The QCD running coupling 

In the QCD lagrangian the only parameters with the dimension of energy are the quark masses 
(in units h = c = 1). Thus at the classical level massless QCD is scale invariant. But the scale 
invariance of massless QCD does not remain true in the quantum theory. The scale symmetry of 
the classical theory is unavoidably destroyed by the regularization and renormalization procedure 
which introduce a dimensional parameter in the quantum version of the theory. When a symmetry 
of the classical theory is necessarily destroyed by quantization, regularization and renormalization 
one talks of an "anomaly". So, in this sense, scale invariance in massless QCD is anomalous. 

While massless QCD is finally not scale invariant, the departures from scaling ai^e asymptot- 
ically small, logarithmic and computable. In massive QCD there are additional mass corrections 
suppressed by powers of m/E, where E is the energy scale (for processes that are non singular in 
the limit m — )• 0). At the parton level (q and g) we can conceive to apply the asymptotic predictions 
of massless QCD to processes and observables (we use the word "processes" for both) with the 
following properties ("hard processes"), (a) All relevant energy variables must be large: 

Ei = ZiQ, Q»mf, z,: scaling variables o(l) (2.1) 

(b) There should be no infrared singularities (one talks of "infrared safe" processes), (c) The pro- 
cesses concerned must be finite for m (no mass singularities). To possibly satisfy these criteria 
processes must be as "inclusive" as possible: one should include all final states with massless gluon 
emission and add all mass degenerate final states (given that quarks are massless also q — q pairs 
can be massless if "coUinear", that is moving together in the same direction at the common speed 
of hght). 

In perturbative QCD one computes inclusive rates for partons (the fields in the lagrangian, that 
is, in QCD, quarks and gluons) and takes them as equal to rates for hadrons. Partons and hadrons 
are considered as two equivalent sets of complete states. This is called "global duality" and it is 
rather safe in the rare instance of a totally inclusive final state. It is less so for distributions, like 
distributions in the invariant mass M ("local duality") where it can be reliable only if smeared over 
a sufficiently wide bin in M. 

For the derivation of the running coupling the basic framework is the renormalization group 
formahsm. In this section we denote with a either a, in the case of QCD or its analogue in QED. 
For a sufficiently well defined Green function G, the renormalization group equation (RGE) in 
massless QCD can be written as: 

+ + 7o(«)]G.„ = (2.2) 

where a is the renormalized coupling defined at the scale /i (jj. cannot be zero because of infrared 
divergences) and 

^ da 
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and 

7o(a) = (2.4) 
5 log 

Note that /3(a) does not depend on which Green function G we are considering; actually it is a 
property of the theory and of the renonnalization scheme adopted, while Yg{(^) also depends on 
G. Strictly speaking the RGE as written above is only valid in the Landau gauge (A = 0). In other 
gauges an additional term that takes the variation of the gauge fixing parameter A should also be 
included. We omit this term, for simplicity, as it is not relevant at the 1-loop level. 

Assume that we want to apply the RGE to some hard process at a large scale Q, related to a 
Green function G that we can always take as dimensionless (by multiplication by a suitable power 
of Q). Since the interesting dependence on Q will be logarithmic we introduce the variable t as : 

? = log^ (2.5) 

Then we can write Gren = F{t, a,Xi) where xi are scaling variables (we often omit to write them in 
the following). In the naive scaling limit F should be independent of t, according to the classical 
intuition that massless QCD is scale invariant. To find the actual dependence on ?, we want to solve 
the RGE 

[-|- + i3(a)^ + 7g(«)]G,,„ = (2.6) 
ot da 

with a given boundary condition at t = (or = pL^): F(0,a). The solution of this general 
equation is given by: 

F{t,a) = F[0,a(0]exp / ^j-r,doc' (2.7) 



where the "running coupling" a{t) is defined by: 

.a{f) I 



-da' (2.8) 



/3(a') 

Note that from this definition it follows that a(0) = a, so that the boundary condition is also 
satisfied. To prove that F(?,a) in eq. 2^ is indeed the solution, we first take derivatives with 
respect of t and a (the two independent variables) of both sides of eq. (2^). By taking d/dt we 
obtain 

1 d0C{t) 

^ ~ I5{a{t) dt ^ ^ ^ 



We then take d/da and obtain 



1 1 da(t) 

j3(a) j8(a(0 da 



These two relations make explicit the dependence of the running coupling on t and a 

da{t) 



dt 

dajt) _ pjajt)) 
da j3(a) 



i3(«(0) (2.11) 

(2.12) 
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Using these two equations one immediately checks that F{t,a) is indeed the solution. In fact the 
sum of the two derivatives acting on the factor F[0, cc{t)] vanishes and the exponential is by itself a 
solution of the complete equation. Note that the boundary condition is also satisfied. The important 
point is the appearance of the running coupling as the quantity that determines the asymptotic 
departures from scaling. 

The next step is to study the functional form of the running coupling. From eq. ( 2.1 1[ ) we see 



that the rate of change with t of the running coupling is determined by the j8 function. In turn j8 (a) 



is determined by the jx dependence of the renormahzed coupling through eq. ( |2.3[ ). Clearly there is 
no dependence on jj. of the basic 3-gluon vertex in lowest order (order e). The dependence starts at 
1-loop, that is at order (one extra gluon has to be emitted and reabsorbed). Thus we obtain that 
in perturbation theory: 



(2.13) 



5 log /1 2 



Recalling that a = /An, we have: 

da 



2g -, , oc / oc (2.14) 



Thus the behaviour of j8(a) in perturbation theory is as follows: 

jS(a) = ±ba?[l + b'a + ...] (2.15) 

Since the sign of the leading term is crucial in the following discussion, we stipulate that always 
b>0 and we make the sign explicit in front. 
By direct calculation at 1-loop one finds: 

QED: j3(a) ~ +ha^ + b = V ^^^^ (2.16) 

^ 371 

where Nc = 3 for quarks and Nc = I for leptons and the sum runs over all fermions of charge Q/e 
that are coupled. Also, one finds: 

QCD: /3(a) ~ -ha^ + b = ^^Nc-2nf ^^jv) 

1271 

where, as usual, «/ is the number of coupled (see below) flavours of quarks (we assume here 
that nf < 16 so that b > in QCD). If a{t) is small we can compute j8(a(f)) in perturbation 
theory. The sign in front of b then decides the slope of the coupling: a{t) increases with t (or 
Q^) if j3 is positive at small a (QED), or a{t) decreases with t (or Q^) if jS is negative at small 
a (QCD). A theory like QCD where the running coupling vanishes asymptotically at large is 



called (ultraviolet) "asymptotically free". An important result that has been proven [ ]12[ ] is that in 4 
spacetime dimensions all and only non-abelian gauge theories ai^e asymptotically free. 

Going back to eq. ( ^ ) we replace j8(a) ~ ±ba^, do the integral and perform a simple 
algebra. We find 

QED: a(t) ~ — ^ — (2.18) 
l-b(Zt 



and 
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A slightly different forni is often used in QCD. Defining Aqcd by 1/a = blogjj. /Aq^jj, we can 
write: 

«(0 ~ -J-^ = T-^ = (2-20) 

i+bt blog + Z;log g blog 

The parameter pL has been traded for the parameter Aqcd- We see that a{t) decreases logarith- 
mically with and that one can introduce a dimensional parameter Aqcd that replaces jLt. Often 
in the following we will simply write A for Aqcd- Note that it is clear that A depends on the 
particular definition of a, not only through the defining scale /i but also on the renormalization 
scheme adopted. It also depends on the number «/ of coupled flavours through the parameter b, 
and in general through the j8 function. It is very important to note that QED and QCD aie theories 
with "decoupling": up to the scale Q only quarks with masses m << Q contribute to the running 
of a. This is clearly very important, given that all applications of perturbative QCD so far apply to 



energies below the top quark mass m,. For the validity of the decoupling theorem [ |13[ ] it is neces- 
sary that the theory where all the heavy particle internal lines are eliminated is still renormalizable 
and that the coupling constants do not vary with the mass. These requirements are true for heavy 
quaifc in QED and QCD, but aie not true in the electroweak theory where the elimination of the 
top would violate SU (2) symmetry (because the t and b left-handed quarks ai^e in a doublet) and the 
quark couplings to the Higgs multiplet (hence to the longitudinal gauge bosons) are proportional 
to the mass. In conclusion, in QED and QCD, quarks with m » Q do not contribute to tif in 
the coefficients of the relevant j3 function. The effects of heavy quarks are power suppressed and 
can be taken separately into account. For example, in e^e^ annihilation for Inic < Q < 2mh the 
relevant asymptotics is for «y = 4, while for 2mb < Q < 2m, rif = 5. Going accross the b threshold 
the p function coefficients change, so the a{t) slope changes. But a{t) is continuous, so that A 
changes so as to keep a{t) constant at the matching point at 2 ~ o{2mi,). The effect on A is large: 
approximately Aj ~ O.65A4 where A4 5 are for «/ = 4,5. 



Note the presence of a pole in eqs.( |2.18y2.19D at ±bat = 1, called the Landau pole, who 



realized its existence in QED already in the '50's. For /i ~ nie, in QED, the pole occurs beyond 
the Planck mass. In QCD the Landau pole is located for negative t or ai Q < pi in the region of 
light hadron masses. Clearly the issue of the definition and the behaviour of the physical coupling 
(which is always finite, when defined in terms of some physical process) in the region around the 
perturbative Landau pole is a problem that lies outside the domain of perturbative QCD. 

The non leading terms in the asymptotic behaviour of the running coupling can in principle 



be evaluated going back to eq. ( 2.15 ) and computing b' at 2-loops and so on. But in general the 
perturbative coefficients of j8(a) depend on the definition of the renormalized coupling a (the 
renormalization scheme), so one wonders whether it is worthwhile to do a complicated calculation 
to get b' if then it must be repeated for a different definition or scheme. In this respect it is interest- 
ing to remark that actually one can easily prove that, due to the definition of j8(a), both b and b' 
are independent of the definition of a, while higher order coefficients do depend on that. 
In QCD {Nc = 3) one has obtained @: 
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By taking b' into account one can write the expression of the running coupling at next to the leading 
order (NLO): 

«(e') = «Lo(e')[l - b'aLo{Q^)\og\og^ + ...] (2.22) 

where a^^ = bXogQ^ / is the LO result (actually at NLO the definition of A is modified ac- 
cording to b log pi^ /A^ = l/a + b'\og ba). 

At present the universally adopted definition of a, is in terms of dimensional regularization 
[p4|l, with the so-called Modified Minimal Subtraction (MS) scheme 1 15] and a value quoted for 



is normally referring to this definition. Different measurements obtained at different energy scales 
are usually made comparable to each other by translating the result in temis of the value at the Z 
mass: as{mz). 

The third ^ and fourth [|18]| coefficients of the QCD jS function are also known in the MS 



prescription (recall that only the first two coefficients are scheme independent). The calculation of 
the last term involved the evaluation of some 50,000 4-loop diagrams. Translated in numbers, for 
rif = 5 one obtains : 

p{a) = -0.610a2[l + 1.261...- + 1.475. ..{-f + 9.836.. (2.23) 

71 71 71 

It is interesting to remark that the expansion coefficients are of order 1 or 10 (only for the last one), 
so that the MS expansion looks reasonably well behaved. 

Summarizing, massless classical QCD is scale invariant. But the procedure of quantization, 
regularization and renormalization necessarily breaks scale invariance. In the quantum QCD theory 
there is a scale of energy. A, which, from experiment, is of the order of a few hundred MeV, 
its precise value depending on the definition, as we shall see in detail. Dimensionless quantities 
depend on the energy scale through the running coupling which is a logarithmic function of Q^/A^. 
In QCD the running coupling decreases logarithmically at large (asymptotic freedom), while in 
QED the coupling has the opposite behaviour. 

3. Measurements of a. 

Very precise and reUable measurements of as{mz) are obtained from e^e^ colliders (in partic- 
ular LEP), from deep inelastic scattering and from the hadron Colliders (Tevatron and LHC). The 



"official" compilation due to Bethke pO| ] and included in the 2012 edition of the PDG [ p^ leads 
to the world average as{mz) = 0.1 184 it 0.0007. A similar analysis with equivalent conclusions 
was recently presented in ref. [0. The agreement among so many different ways of measuring 
Us is a strong quantitative test of QCD. However for some entries the stated error is taken directly 
from the original works and is not transparent enough as seen from outside (e.g. the lattice deter- 
mination). In my opinion one should select few theoretically simplest processes for measuring 
and consider all other ways as tests of the theory. Note that this is what is usually done in QED 
where a is measured from one single very precise and theoretically reliable observable (one pos- 



sible calibration process is at present the electron g-2 [22]). The cleanest processes for measuring 
Us are the totally inclusive ones (no hadronic coixections) with light cone dominance, like Z decay, 
scaling violations in DIS and perhaps T decay (but, for T, the energy scale is dangerously low). We 
will review these cleanest methods for measuring a, in the following. 
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3.1 ttv from e^e colliders 

The totally inclusive processes for measuring a, at e^e^ colliders are hadronic Z or T decays 
(/?/ T = r(Z, T — )• hadrons) /T{Z, T — leptons), Oh, Oj, [Op = 127rr/rf / {m\T\) with F=h or 1 are 
the hadronic or leptonic cross sections at the Z peak], and the total Z width Yz = SF/ + F/, + F,„v. 
For each of this quantities, for example 7?/, one can write a general expression of the form: 

R, ~ R^'^il + 5qcd + 8np) (3.1) 

where 7?^^ is the electroweak-corrected Bom approximation, 5qcd, ^np are the perturbative (log- 
arithmic) and non perturbative (power suppressed) QCD corrections. For a measurement of as (we 
always refer to the MS definition of a^) at the Z resonance peak one can use all the information 
from Ri, Tz and a/,,;. In the past the measurement from 7?/ was preferred (by itself it leads to 
Oli(wi2) = 0. 1226 ± 0.0038, a bit on the large side) but after LEP there is no reason for this prefer- 
ence. In all these quantities enters through F/,, but the measurements of, say, Tz, Ri and a/ are 
really independent as they are affected by an entirely different systematics: Fz is extracted from the 
line shape, /?; and (7; are measured at the peak but Ri does not depend on the absolute luminosity 
while Oi does. The most sensitive single quantity is a/. It gives as{m.z) = 0.1 183 ±0.0030. The 
combined value from the measurements at the Z (assuming the validity of the SM and the observed 



Higgs mass) is []23|]: 

a,(mz) =0.1187 ±0.0027 (3.2) 
By adding all other electroweak precision electroweak tests (in particular m-w) one similarly finds 



[24] 



a, („jz) =0.1 186 ±0.0026 (3.3) 

These results have been obtained from the Sqcd expansion up to and including the C3 term of order 
. But by now the C4 term (NNNLO!) has also been computed [|^] for inclusive hadronic Z and 
T decay. This remarkable calculation of about 20.000 diagrams, for the inclusive hadronic Z width, 
led to the result, for = 5 and a, = as{mz) / 7i: 

8qcd = [1 ± a. ± 0.76264 a] - 15.49 a] - 68.2 a] + ...] (3.4) 



This result can be used to improve the value of as{niz) from the EW fit given in eq. (3.3) that 
becomes: 

a,(mz) =0.1 190 ±0.0026 (3.5) 

Note that the error shown is dominated by the experimental errors. Ambiguities from higher pertur- 
bative orders [^], from power corrections and also from uncertainties on the Bhabha luminometer 
(which affect a/,,/) p7| ] are very small. In particular, having now fixed ni// does not decrease the 



error significantly [ [28[ , [29[ ]. The main source of error is the assumption of no new physics, for 
example in the Zbb vertex that could affect the F/, prediction. 

We now consider the measurement of as{mz) from T decay. R^ has a number of advantages 
that, at least in part, tend to compensate for the smallness of = 1 .777 GeV. First, 7?^ is maximally 
inclusive, more than Rg+e- {s), because one also integrates over all values of the invariant hadronic 
squared mass: 

1 ds , s 



R, = - -^{l-^yimn,{s) (3.6) 
n Jo mj. 
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As we have seen, the perturbative contribution is now known at NNNLO [25]. Analyticity can be 
used to transform the integral into one on the circle at \s\ = m^: 

R^ = ^-i ,^(i-4)'n-(^) (3.7) 

2ni J\s\=m\ m% m\ 

Also, the factor (1 — is important to kill the sensitivity the region Re\s\ = where the phys- 
ical cut and the associated thresholds are located. Still the sensitivity to hadronic effects in the 
vicinity of the cut is a non negligible source of theoretical eiTor that the formulation of duality 
violation models try to decrease. But the main feature that has attracted attention to x decays for 
the measurement of as{niz) is that even a rough determination of Aqcd at a low scale Q ~ 
leads to a very precise prediction of at the scale niz, just because in logQ/Agcz) the value of 
Aqcd counts less and less as Q increases. The absolute error on shrinks by a factor of about 
one order of magnitude going from o;.v(mT) to as{mz). Still I find a httle suspicious that to obtain 
a better measurement of as{rnz) you have to go down to lower and lower energy scales. And in 
fact, in general, in similar cases one finds that the decreased control of higher order perturbative 
and of non perturbative corrections makes the apparent advantage totally illusory. For a, from 7?^ 
the quoted amazing precision is obtained by taking for granted that coixections suppressed by 1/m^ 
are negligible. The argument is that in the massless theory, the light cone expansion is given by: 

ZERO <0,> <0e> 

ONP = J h C4 7 h C6 7 \ (3.8) 

jrq. in^ 

In fact there are no dim-2 Lorentz and gauge invariant operators. For example, Y^AsfiS^^^ where 
are the gluon fields, is not gauge invariant. In the massive theory, the ZERO is replaced by 
the light quark mass-squared m^. This is still negligible if m is taken as a lagrangian mass of a 
few MeV. If on the other hand the mass were taken to be the constituent mass of order Aqcd, this 
term would not be negligible at all and would substantially affect the result (note that as{ni^)/n ~ 
0.1 ~ (0.6 GeV/m^f and that Aqcd for 3 fl avours is large). The principle that coefficients in 
the operator expansion can be computed from the perturbative theory in terms of parton masses has 
never been really tested (due to ambiguities on the determination of condensates) and this particular 
case where a ZERO appears in the massless theory is unique in making the issue crucial. Many 
distinguished people believe the optimistic version. I am not convinced that the gap is not filled 
up by ambiguities of 0(Ag(^^/wi^) from 5pen- the [ZERO/m^] terms in eq. 3.8 are vulnerable 



to possible ambiguities from the so-called ultraviolet renormalons [30|. In fact one must keep in 
mind that the QED and QCD perturbative series, after renormalization, have all their coefficients 
finite, but the expansion does not converge. Actually the perturbative series is not even Borel 
summable (for reviews, see, for example refs. [^). After the Borel resummation, for a given 
process one is left with a result which is ambiguous by terms typically down by exp—n/{ba), 
with n an integer and b is the absolute value of the first p function coefficient. In QED these 
corrective terms are extremely small and not very important in practice. On the contrary in QCD 
o; = CCs{Q^) ~ l/{b\ogQ^ /A^) and the ambiguous teniis ai^e of order (1/2^)", that is are power 
suppressed. It is interesting that, through this mechanism, the perturbative version of the theory is 
able to somehow take into account the power suppressed coixections. A sequence of diagrams with 
factorial growth at large order n is made up by dressing gluon propagators by any number of quark 
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bubbles together with their gauge completions (renormalons).The problem of the precise relation 
between the ambiguities of the perturbative expansion and the higher twist corrections has been 
discussed in recent years [|T]]. 

There is a vast and sophisticated literature on from T decay. Unbelievably small eixors on 
CCs{niz) ai"e obtained in one or the other of several different procedures and assumptions that have 
been adopted to end up with a specified result. With time there has been an increasing awareness 
on the problem of controlling higher orders and non perturbative effects. In particular fixed order 
perturbation theoiy (FOPT) has been compared to resummation of leading beta function effects 
in the so called contour improved perturbation theory (CIPT). The results are sizably different in 
the two cases and there have been many arguments in the literature on which method is best. One 
important progress comes from the experimental measurement of moments of the T decay mass 



distributions, defined by modifying the weight function in the integral in eq.(|3.q). In principle 
one can measure from the sum rules obtained from different weight functions that emphasize 
different mass intervals and different operator dimensions in the light cone operator expansion. A 
thorough study of the dependence of the measured value of on the choice of the weight function 



and in general of higher order and non perturbative corrections has appeared in ref . [|32[] and I advise 
the interested reader to look at that paper and the references therein. It would be great to be able to 
fit both the value of a, and of the coefficient of the term at the same time. Unfortunately a 
good precision on is only obtained if that coefficient is fixed. 

We consider here the recent evaluations of from T decay based on the NNNLO perturbative 



calculations [25] and different procedures for the estimate of all sorts of coixections. From the 



papers given in refs. [33] we obtain an average value and eixor that agrees with the Erler and 
Langacker values given in the PDG'12 [[l9|]: 

a,{m^) = 0.3285 ±0.018 (3.9) 

or 

a, („iz) =0.1 194 ±0.0021 (3.10) 

In any case, one can discuss the error, but what is true and remarkable, is that the central value of 
from T decay, obtained at very small Q^, is in good agreement with all other precise determinations 
of as at more typical LEP values of Q 



2 



3.2 as from Deep Inelastic Scattering 

In principle DIS is expected to be an ideal laboratory for the determination of as but in practice 
the outcome is still to some extent open. QCD predicts the dependence of F{x,Q^) at each 
fixed X, not the x shape. But the dependence is related to the x shape by the QCD evolution 
equations. For each x-bin the data allow to extract the slope of an approximately straight line in 
dlogF {x,Q^) / dlogQ^: the log slope. The span and the precision of the data are not much 
sensitive to the curvature, for most x values. A single value of Agco must be fitted to reproduce 
the collection of the log slopes. For the determination of a, the scaling violations of non-singlet 
structure functions would be ideal, because of the minimal impact of the choice of input parton 
densities. We can write the non-singlet evolution equations in the form: 

|/o,F(x,,) = ^/';^^^?„(^,a,«)) (3^11, 

at 271 Jx y r(x,t) y 



10 



Measurements of 



Guido Altarelli 



where Pqq is the splitting function. At present NLO and NNLO corrections are known. It is cleai^ 
from this fomi that, for example, the normalization error on the input density drops away, and the 
dependence on the input is reduced to a minimum (indeed, only a single density appears here, while 
in general there are quai^k and gluon densities). Unfortunately the data on non-singlet structure 
functions are not very accurate. If we take the difference of data on protons and neutrons, Fp — Fn, 
experimental errors add up in the difference and finally are large. The F-^vn data are directly non- 
singlet but are not very precise. Another possibility is to neglect sea and glue in F2 at sufficiently 
large x. But by only taking data &tx> xq one decreases the sample, introduces a dependence on xq 
and an error from residual singlet terms. A recent fit to non singlet structure functions in electro- or 
muon-production extracted from proton and deuterium data, neglecting sea and gluons ntx> 0.3 



(error to be evaluated) has led to the results [ 34 1 



as{mz) = 0.lU%±0mi9{exp) + 1 {NLO) (3.12) 
a,(mz) = 0.1134±0.0020(e;cp) + ? {NNLO) (3.13) 

The central values are rather low and there is not much difference between NLO and NNLO. The 
question marks refer to the uncertainties from the residual singlet component at ;c > 0.3 and also 
to the fact that the old BCDMS data, whose systematics has been questioned, are very important at 
X > 0.3 and push the fit towards small values of a^. 

When one measures from scaling violations in F2 measured with e or beams, the data 
are abundant, the statistical eiTors are small, the ambiguities from the treatment of heavy quarks 
and the effects of the longitudinal structure function Fl can be controlled, but there is an increased 
dependence on input parton densities and especially a strong correlation between the result on 
and the adopted parametrization of the gluon density. In the following we restrict our attention to 
recent determinations of a, from scaling violations at NNLO accuracy, as, for example, those in 



refs. [|35|, ^ that report the results, in the order: 



as{mz) = 0.1 134 ± 0.001 l(£';c/7) + ? (3.14) 
a,(mz) = 0.1 158 ±0.0035 (3.15) 

In the first line my question mark refers to the issue of the ttv-gluon conelation. In fact tends to 
slide towards low values (a^ ~ 0.113 — 0.116) if the gluon input problem is not fixed. Indeed, in 
the second line from ref. [^, the large error also includes an estimate of the ambiguity from the 
gluon density parametrization. One way to restrict the gluon density is to use the Tevatron high pr 
jet data to fix the gluon parton density at large x that, via the momentum conservation sum rule, 
also constrain the small x values of the same density. Of course in this way one has to go outside 
the pure domain of DIS. Also, the jet rates have been computed at NLO only. In a simultaneous fit 
of as and the parton densities from a set of data that, although dominated by DIS data, also contains 



Tevatron jets and Drell- Yan production, the result was [37| 



a^(mz) =0.1171 ±0.0014 + ? (3.16) 



The authors of ref. [37] attribute their larger value of to a more flexible parametrization of 



the gluon and the inclusion of Tevatron jet data that are important to fix the gluon at large x. An 
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alternative way to cope with the gluon problem is to drastically suppress the gluon parametrization 



rigidity by adopting the neural network approach. With this method, in ref. [|38|], from DIS data 
only, treated at NNLO accuracy, the following value was obtained: 

a,(mz) =0.1 166 ±0.0008(e;cp)±0.0009(f/t) +? (3.17) 

where the stated theoretical error is that quoted by the authors within their framework, while the 
question mark has to do with possible additional systematics from the method adopted. Interest- 
ingly, in the same approach, by also including the Tevatron jets and the Drell-Yan data not much 
difference is found: 



a,{mz)=OAm±0.mQl{exp)±0.mm{th) +? (3.18) 

We see that when the gluon input problem is suitably addressed the fitted value of is increased. 

As we have seen there is some spread of results, even among the most recent determinations 
based on NNLO splitting functions. We tend to favour determinations from the whole DIS set of 
data (i.e. beyond the pure non singlet case) and with attention paid to the gluon ambiguity problem 
(even if some non DIS data from Tevatron jets at NLO have to be included). A conservative pro- 
posal for the resulting value of Us from DIS, that emerges from the above discussion is something 
like: 

a, (mz) =0.1165 ±0.0020 (3.19) 

The central value is below those obtained from Z and T decays but perfectly compatible with those 
results. 



3.3 Other av(mz) Measurements as QCD Tests 

There are a number of other determinations of that are important because they arise from 
qualitatively different obsei^vables and methods. All together they provide an impressive set of 
QCD tests. Here I will give a few examples of these interesting measurements. 

A classic set of measurements is from a number of infrared safe observables related to event 
rates and jet shapes in e^e^ annihilation. One important feature of these measurements is that they 
can be repeated at different energies in the same detector, like the JADE detector in the energy 
range of PETRA or the LEP detectors from LEPl to LEP2 energies. As a result one obtains a 
striking direct confirmation of the running of the coupling according to the renormalization group 
prediction. The perturbative part is known at NNLO [|9|| and resummations of leading logs arising 
from the vicinity of cuts and/or boundaiies have been performed in many cases using effective 
field theory methods. The main problem of these measurements is the possible large impact of non 
perturbative hadronization effects on the result and therefore on the theoretical eixor According 
to ref.[pO|] a summary result that takes into account the central values and the spread from the 
JADE measurements, in the range 14 to 46 GeV, at PETRA is given by: as{mz) = 0.1172 it 
0.0051, while from the ALEPH data at LEP, in the range 90 to 206 GeV, the reported value Q 
is a^(mz) = 0.1224 lb 0.0039. It is amazing to note that among the related works there ai^e a 
couple of papers by Abbate et al [^, 12 1 where an extremely sophisticated formalism is developed 
for the thrust distribution, based on NNLO perturbation theory with resummations at NNNLL 
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plus a data/theoiy based estimate of non peiturbative corrections. The final quoted results ai^e 
unbelievably precise: as{mz) = 0.1135 ±0.0011 from the tail of the Thrust distribution pT] ] and 



0Csi>nz) = 0.1140 ± 0.0015 from the first moment of the Thrust distribution [ p2| | (note the low 
central value). I think that this is a good example of an underestimated error which is obtained 
within a given machinery without considering the limits of the method itself. Another allegedly 



very precise determination of as{mz) is obtained from lattice QCD by several groups [Q with 
different methods and compatible results. A value that summarizes these different results is [19] 
o;.v('Mz) = 0.1 185 ± 0.0007. With all due respect to lattice people I think this small eixor is totally 
umplausible. 

4. Conclusion: My Recommended Value of as{mz) 

According to my proposal to calibrate (Xs{mz) from the theoretically cleanest and most trans- 
parent methods, identified as the totally inclusive, light cone operator expansion dominated pro- 
cesses, I collect here my understanding of the results: from Z decays and EW precision tests. 



eq.(y): 

a^.(mz) =0.1 190 ±0.0026; (4.1) 
from scahng violations in DIS, eq.( |3.19| ): 

a,{mz) = 0.1165 ± 0.0020; (4.2) 



from /?i:,eq.(3.10) 



a,(mz) =0.1 194 ±0.0021. (4.3) 
If one wants to be on the safest side one can take the average of Z decay and DIS : 

a,(mz) =0.1 174 ±0.0016. (4.4) 

This is my recommended value. If one adds to the average the relatively conservative 7?^ value and 
error given above in eq. 3.10| , that takes into account the dangerous low energy scale of the process. 



one obtains: 

a^(mz) = 0.1 184 ±0.0011. (4.5) 
Note that this is essentially coincident with the "official" average with a moderate increase of the 



error Thus we see that a sufficiently precise measure of aj(mz) can be obtained, eqs. (1.4,4.5), 
by only using the simplest processes where the control of theoretical errors is maximal. One is left 
free to judge whether a further restriction of theoretical errors is really on solid ground. 
The value of A (for Uf = 5) which coiTesponds to eq. ( p] ) is: 

A5 =202±18MeV (4.6) 



while the value from eq. ( [4.5| ) is: 

As =213±13MeV (4.7) 
In conclusion, we see that a reasonably precise measure of a^(mz) can be obtained, eqs. 



( p|J4.5[ ), by only using the simplest processes where the control of theoretical enws is maxi- 
mal. One is left free to judge whether a further decrease of the eiTor is really worthwhile at the 
price of a loss in transparency and rigour 



13 



Measurements of 



Guido Altarelli 



References 

[1] M. Gell-Mann, Acta Phys. Austriaca Suppl. IX (1972) 733; H. Fritzsch and M. Gell-Mann, Proc. XVI 
Int. Conf. on High Energy Physics, Chicago-Batavia, 1972; H. Fritzsch, M. Gell-Mann and H. 
Leutwyler, Phys. Lett. B 47, 365 (1973) 

[2] D. Gross and F Wilczek, Phys. Rev. Lett. 30, 1343 (1973), Phys. Rev. D 8, 3633 (1973); H.D. 
PoUtzer, Phys. Rev. Lett. 30, 1346 (1973) 

[3] S. Weinberg, Phys. Rev. Lett. 31, 494 (1973) 

[4] Y.L. Dokshitzer and V.A. Khoze, Basics of Perturbative QCD, Ed. Frontieres (1991); T. Muta, 
Foundation of Quantum Chromodynamics: an Introduction to Perturbative Methods in Gauge 
Theories, World Sci., 3rd ed. (2009); F. J. Yndurain, The Theory of Quark and Gluon Interactions, 
Springer, 4th ed., (2006). G. Altarelli, The Development of Perturbative QCD, World Scientific, 
(1994); W. Greiner, S. Schramm and E. Stein, Quantum Chromodynamics, Springer, 3rd ed. (2007); 
Handbuch of QCD, ed. by M. Shifman, Vol.1-4, World Sci. (2001); R. K. Ellis, W. J. Stirling, and B. 
R. Webber, QCD and Collider Physics, Cambridge Monographs, (2003); G. Dissertori, I. Knowles 
and M. Schmelling, Quantum Chromodynamics: High Energy Experiments and Theory, Oxford Univ. 
Press, 2003; S. Narison,(2CZ) as a Theory of Hadrons: From Partons to Confinement, Cambridge 
Monographs, (2007); J. C. Collins, Foundations of Perturbative QCD , Cambridge Monographs, 
(2011) 

[5] G. Altarelli, ArXiv: 1303.2842 
[6] A. S. Kronfeld, ArXiv: 1203. 1204 

[7] S. Weinberg, Physica A 96 (1979) 327, J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142; 
Nucl.Phys.B250(1985)465;J. Bijnens, G. Ecker and J. Gasser, hep-ph/941 1232, in The Second 
DA4>NE Physics Handbook (Frascati, 1995); J. Bijnens and U. Meissner, hep-ph/9901381; G. Ecker, 
hep-ph/9805500, hep-ph/001 1026; H. Leutwyler, hep-ph/9609465, hep-ph/0008124; A. Pich, 
hep-ph/9806303; E. de Rafael, hep-ph/9502254; L. Maiani, G. Pancheri and N. Paver, The Second 
DA4>NE Physics Handbook (Frascati, 1995); G. Ecker, Prog. Part. Nucl. Phys. 35 (1995) 1; G. 
Colangelol andG. Isidori, ArXiv :hep-ph/0101264 

[8] J. Gasser, Nucl.Phys.Proc.Suppl.86:257 (2000), ArXiv :hep-ph/99 12548 

[9] E. Eichten and B. Hill, Phys. Lett. B 234, 511 (1990); H. Georgi, Phys. Lett. B 240, 447 (1990); B. 
Grinstein, Nucl. Phys. B 339, 253 (1990); T. Mannel, W. Roberts, and Z. Ryzak, Nucl. Phys. B 368, 
204(1992) 

[10] M. Neubert, Phys. Rept. 245, 259 (1994) ArXiv;hep-ph/9306320; A.V. Manohar and M.B. Wise, 
Camb. Monogr Part. Phys. Nucl. Phys. Cosmol. 10, 1 (2000) 

[11] C. W. Bauer, S. Fleming, and M.E. Luke, Phys. Rev. D 63, 014006 (2000), ArXiv:hep-ph/0005275; 
C.W. Bauer et al., Phys. Rev. D 63, 1 14020 (2001) ArXiv:hep-ph/001 1336]; C.W. Bauer and I.W. 
Stewart, Phys. Lett. B 516, 134 (2001) ArXiv:hep-ph/0107001; C.W. Bauer, D. Pirjol, and I.W. 
Stewart, Phys. Rev. D 65, 054022 (2002), ArXiv:hep-ph/0 109045 

[12] S. Coleman and D. Gross, Phys. Rev. Lett. 31 (1973) 851 

[13] T. Applequist and J. Carazzone, Phys. Rev. Dl 1 (1975)2856 

[14] G. 't Hooft and M. Veltman, Nucl. Phys. B44 (1972) 189; C.G. Bollini and J.J. Giambiagi, Nuovo 
Cimento 12B (1972) 20; J.F Ashmore, Nuovo Cimento Lett. 4 (1972) 289; G.M. Cicuta and E. 
Montaldi, Nuovo Cimento Lett 4 (1972) 329 

[15] G. 't Hooft, Nucl. Phys. B61 (1973) 455; W. A. Bardeen, A. J. Buras, D. W. Duke and T Muta, Phys. 
Rev. D18(1978) 3998 



14 



Measurements of 



Guido Altarelli 



[16] W. Caswell, Phys. Rev. Letters 33 (1974) 244; D. R. T. Jones, Nucl. Phys. B75 (1974) 531 
[17] O.V. Tarasov, A.A. Vladimirov and A. Yu. Zharkov, Phys. Lett. 93B(1980)429 

[18] T. van Ritbergen, J. A.M. Vermaseren and S.A. Larin, Phys.Lett. B400(1997)379; see also M. Czakon, 

Nucl. Phys. B710 (2005) 485 
[19] Particle Data Group, J. Beringer et al, Phys. Rev. D86(2012) 010001 
[20] S. Bethke, ArXiv;0908.1135 
[21] A. Pich, ArXiv:1303.2262 

[22] D. Hanneke, S. Fogwell, and G. Gabrielse, Phys. Rev. Lett. 100, 120801 (2008); D. Hanneke, S. 
Fogwell Hoogerheide, and G. Gabrielse, Phys. Rev. A 83, 052122 (201 1) 

[23] S. Kluth, Rept. Prog. Phys. 69 (2006) 1771 

[24] The LEP Electroweak Working Group, http://lepewwg.web.cern.ch/LEPEWWG/ 

[25] PA. Baikov, K.G. Chetyrkin and J.H. Kuhn, Phys. Rev. Lett. 101 (2008) 012002, ArXiv:0801.1821; 
Phys. Rev. Lett. 104 (2010) 132004, ArXiv;1001.3606P; A. Baikov, K. G. Chetyrkin, J. H. Kuhn and 
J. Rittinger, ArXiv:1210.3594 

[26] H. Stenzel, JHEP 0507:0132,2005 

[27] W. de Boer and C. Sander, Phys.Lett.B585:276,2004 

[28] M. Baak et al. Glitter group, ArXiv: 1209.2716 

[29] M. Grunewald, for the LEP EW Group, private communication 

[30] G Altarelli, P Nason and G. Ridolfi, Z.Phys.C68:257,1995 

[31] G. Altarelli, Proceedings of the E. Majorana Summer School, Erice, 1995, Plenum Press, ed. by A. 
Zichichi; M. Beneke and V. M. Braun, in Handbuch ofQCD. ed. by M. Shifman, Vol.3, pag. 1719, 
World Sci. (2001) 

[32] M. Beneke, D. Boito and M. Jamin, ArXiv: 1210.8038 

[33] M. Beneke and M. Jamin, JHEP 0809, 044 (2008), ArXiv:0806.3156; M. Davier et al., Eur. Phys. J. 
C56, 305 (2008), ArXiv:0803.0979; K. Maltman and T. Yavin, Phys. Rev. D78, 094020 (2008), 
ArXiv:0807.0650; S. Narison, Phys. Lett. B673, 30 (2009), ArXiv:0901.3823; I. Caprini and J. 
Fischer, Eur Phys. J. C64, 35 (2009); ArXiv:0906.521 1; S. Menke, ArXiv:0904.1796; A. Pich, 
;ArXiv: 1107. 1123; B.A. Magradze, ArXiv: 11 12.5958; G. Abbas et al., ArXiv: 1202.2672; D. Boito et 
al., Phys. Rev. D84, 1 13006 (201 1), ArXiv: 1 1 10.1 127; D. Boito et al., ArXiv: 1203.3146 

[34] J. Blumlein, H. Bottcher and A. Guffanti, Nucl.Phys.B774: 1 82,2007, ArXiv:0607200 [hep-ph] 

[35] S. Alekhin, J. Blumlein and S. Moch, ArXiv: 1202.2281 [hep-ph] 

[36] P Jimenez-Delgado and E. Reya, Phys. Rev. D 79 (2009) 074023, ArXiv:0810.4274 

[37] A. D. Martin, W. J. Stirling, R. S. Thome and G. Watt, Eur. Phys. J. C 64 (2009) 653, 
ArXiv:0905.3531 

[38] R. D. Ball, et al. Phys. Lett. B 707 (2012) 66, ArXiv: 1 1 10.2483 

[39] A. Gehrmann-De I^idder, T. Gehrmann, E.W.N. Glover and G. Heinrich, Phys. Rev. Lett. 99 (2007) 
132002, ArXiv:0707.1285; JHEP 071 1 (2007) 058, ArXiv:07 10.0346; Phys. Rev. Lett. 100 (2008) 
172001, ArXiv:0802.0813, JHEP 0712 (2007) 094 [ArXiv:07 11.4711 

[40] G. Dissertori et al, JHEP 0908:036,2009, ArXiv:0906.3436 

[41] R. Abbate, M. Fickinger, A. H. Hoang, V. Mateu and 1. W. Stewart, Phys.Rev. D83 (201 1) 074021, 
ArXiv: 1006.3080 

[42] R. Abbate, M. Fickinger, A. H. Hoang, V. Mateu and I. W. Stewart, Phys.Rev. D86 (2012) 094002, 
ArXiv: 1204.5746 



15 



Measurements of 



Guido Altarelli 



[43] C. McNeile et al., [HPQCD Collab.], Phys. Rev. D82, 034512 (2010) ArXiv: 1004.4285; C.T.H. 
Davies et al., [HPQCD Collab., UKQCD Collab., and MILC Collab.], Phys. Rev. Lett. 92, 022001 

(2004) , ArXiv:0304004; Q. Mason et al, [HPQCD and UKQCD Collaborations], Phys. Rev. Lett. 95 

(2005) 052002, hep-lat/0503005; K. Maltman, et al., Phys. Rev. D78, 1 14504 (2008), 
ArXiv:0807.2020; S. Aoki et al., [PACS-CS Collab.], JHEP 0910, 053 (2009), ArXiv:0906.3906; E. 
Shintani et al., [JLQCD Collab.], Phys. Rev. D82, 074505 (2010), ArXiv:1002.0371; B. Blossier et 
al., [ETM Collab.], ArXiv: 1201. 5770 



16 



